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SUMMARY 


Estimation  of  target  position  using  planar  arrays  is 
discussed.  For  high  signal-to-noise  ratio,  an  analytic 
expression  for  the  rms  error  in  position  is  derived  as  a 
function  of  the  array  properties.  This  is  evaluated  for 
various  array  configurations  and  it  is  shown  to  be  approximately 
constant  across  the  plane  of  the  array.  Comparisons  are  made 
with  arrays  in  which  only  pairs  of  elements  are  used  to  estimate 
position,  and  it  is  shown  that  this  may  lead  to  a  degradation  in 
performance. 
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1  INTRODUCTION 

Methods  of  estimating  target  position  and  maintaining  a  track  on  a  target 
have  been  made  conventionally  using  amplitude  (or  phase)  comparison  monopulse 
in  Which  there  are  two  receiver  channels  for  both  azimuth  and  elevation.  This 
method  fails  when  there  is  more  than  one  target  within  the  scene  or  when  addi¬ 
tional  signals  (multipath,  jaaning)  are  present  and  results  in  inaccuracies  in 
tracking.  It  may  be  possible  to  overcome  these  failures  by  using  an  array  of 
receivers,  which  will  allow  the  tracking  of  many  targets  simultaneously  (Webb, 
Hearn  and  Sleigh,  1984).  However,  it  is  the  purpose  of  this  research  note  to 
consider  the  single  target  case  and  to  evaluate  the  tracking  error  as  a  func¬ 
tion  of  target  position  in  the  scene,  array  size  and  signal-to-noise  ratio. 

The  reason  for  the  emphasis  on  the  single  target  situation  is  that  it  is  likely 
to  occur  in  many  practical  radar  situations  and  it  is  important  to  understand 
fully  the  properties  of  an  array  in  tracking  single  targets.  The  multi-target 
situation  has  been  considered  for  a  linear  array  only  in  RSRE  Memorandum  3587 
and  results  for  a  planar  array  will  be  the  subject  of  a  separate  report.  In 
the  following  section  the  equations  governing  the  problem  are  described  and  a 
pair  of  coupled  equations  for  the  azimuth  and  elevation  positions  of  the  target 
are  derived.  Section  4  derives  an  analytic  expression  for  the  root  mean  square 
tracking  errors  and  this  is  used  in  section  6,  together  with  results  of  a  simu¬ 
lation,  to  illustrate  the  tracking  ability  of  various  array  configurations. 
Finally,  we  conclude  with  suggestions  for  further  work  and  an  example  of  multi¬ 
target  processing. 

The  main  results  of  this  memorandum  show  that  whilst  the  rms  error  in 
target  position  is  not  reduced  significantly  by  employing  an  array  of  receivers 
as  opposed  to  pairs  of  receivers,  the  advantages  of  an  array  are  in  allowing 
accurate  target  tracking  over  a  wider  field  of  view  and  in  the  ability  to  track 
many  targets  simultaneously.  Also,  arrays  based  on  hexagonal  structures  offer 
comparable  performance  to  those  based  on  square  lattices  but  with  a  reduced 
density  of  elements. 


2  ANALYSIS 


Consider  an  array  of  N  receivers,  each  frith  point-spread  function  as  a 
function  of  azimuth  and  elevation  P£(0,$)  assumed  real.  Denote  the  measured 
output  of  the  array  by  Y  ■  (y^  . . .  yN)  where  y£  is  the  complex  output  of 
receiver  i.  For  a  seen?  consisting  of  M  targets  of  complex  amplitudes  Aq 
(a  -  1,  ...M)  at  positions  -  (6  ,$m),  the  error,  e,  between  the  measured 
output,  Y,  and  the  outputs  for  such  a  scene  is  given  by 


e. 

l 


H 

ZA  p. (x  )  -  y. 
m  l  — m  'i 

m*i 


i  -  1,  ...  N  (1) 


For  a  single  target,  M*  1,  of  amplitude  A*Aj.+  iA£,  we  seek  a  solution  which 
minimises  the  square  error,  e  , 


e2  *  «V  *  ^  CCyJ  -  Arp.(6,^))2+  (yj  -  A.p.«M))2}  (2) 

i*l 

where  y|,  y^  are  the  real  and  imaginary  parts  of  y£  (in  phase  and  quadrature) 
and  pj(§,$)  has  been  assumed  real. 

Differentiating  e?  with  respect  to  the  unknowns  Aj,  A£,  6  and  4>,  and 
equating  to  zero  gives  a  set  of  four  equations  in  four  unknowns.  Eliminating 
Ay  and  A£  between  these,  results  in  a  pair  of  coupled  equations  for  6  and  <t> 


and 


4?  piy')2  4(^  piyj)2)] 

p4?p42  *CSp4H 


(3) 


-  0 


(4) 


The  solution  to  these  equations  depends  on  the  number  of  receiver  elements,  N 
their  point-spread  functions,  p,  and  the  signal-to-noise  ratio.  In  general, 
it  will  be  necessary  to  solve  these  equations  numerically.  However,  for  high 
signal-to-noise  ratio,  it  is  possible  to  obtain  analytic  expressions  for  the 
rms  tracking  errors  in  the  6  and  $  directions.  Firstly,  we  give  a  definition 
of  signal-to-noise  ratio  for  an  array. 
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SIGNAL-TO-NOISE  RATIO  DEFINITION 


It  is  not  inmediately  obvious  what  the  definition  of  signal-to-noise  ratio 
should  be  for  an  array  of  receivers.  For  a  given  scene,  the  signal  power,  S, 
could  be  taken  as 


where  Y  is  the  image  vector,  and  the  noise  power,  N,  as 
N 

n  -  £J<|yil2> 

with  measured  when  there  are  no  targets  present  in  the  scene. 

The  difficulty  with  the  above  definitions  arises  when  comparison  of  results 
relating  to  different  array  sizes  is  made  on  the  basis  of  equal  signal-to-noise 
ratio  since  both  definitions  are  themselves  functions  of  array  size. 

It  is  proposed,  therefore,  that  comparisons  be  made  between  arrays  of 
different  si2e  using  a  given  signal-to-noise  which  depends  on  the  scene  and  the 
point-spread  functions  of  the  elements  only.  For  a  single  target  in  the  scene, 
this  is  defined  as  the  maximum  value  of  signal-to-noise  ratio  on  a  single 
element  as  the  target  is  moved  across  all  points  in  the  scene.  This  definition 
enables  sub-arrays  of  the  total  array  to  be  processed  without  having  to  redefine 
the  signal-to-noise  ratio  each  time. 

4  ANALYTIC  SOLUTIONS 


Suppose  that  the  position  of  the  single  source  is  at  (60,4>0)  and  the  image 
is  corrupted  by  noise,  n,  so  that 


Y  -  I  +  n  , 


(5) 


where  is  the  noiseless  image  given  by 

I  .  »  p.  (6  ,  ♦  )  (6) 

ox  *i'  o  ’  o 

and  n  is  the  noise  vector  added  to  the  image  satisfying 

<|nj|2>  -  < |nT | 2>  -  o2/2  (7) 

For  a  point-spread  function  having  its  maximum  at  (0,0)  the  signal-to-noise 
ratio  is 


p2(0,0)/o2 

For  o2/p2  (0,0)  «  1,  we  look  for  solutions  of  (3)  and  (4)  for  (8,$)  close  to 
(0O»*O**  Expanding  pj(e,$)  in  powers  of  ee  -  (0-0o),  *  ($-*0)  using  Taylor's 

theorem  gives 


t-r  * 

L  * 
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*2fj 

+  2eee4>  3634) 


Using  equations  (5),  (6)  and  (8)  and  substituting  into  equations  (3)  and  (4) 
gives  a  pair  of  simultaneous  equations  for  £q  and 


*lce  *  bi%  *  ci  *  0 

Ve  *  b2%  *  c2  ■  0 


where 


1  Vv  3  3  V*. 

1  6o’V  3  3  V>o 

J  0  0  J 

(?w*M  tJ } 


Similarly  for  a2,  b2  and  <c2>. 

Solving  for  eQ  and  (ajb2  -  a2bj  #  0)  gives: 
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IfT  “  Bh. } 


B(he*  “  bee11^^ +  feh*$  “  2f*f  ehe^  +  f$h 


with  <e^>t  defined  as  <Eq>  ♦  <e|>,  given  by 


(f!  +  f2-B(haa+haJ) 


B(he$  “  heeh*<r  +  f eh*4>  ~  2f$f ehe$  +  f4» 
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The  expressions  for  the  rms  errors  given  here  are  valid  for  any  array  configura¬ 
tion  and  any  element  point-spread  function  in  the  limit  of  large  signal-to-noise 
ratio.  For  an  array  with  elements  symmetrically  placed  in  the  6  and  <t>  directions 
about  the  origin  and  each  element  with  a  symmetrical  point-spread  function, 
equations  (10)  and  (11)  simplify  considerably  for  a  target  on  boresight  to  give 


<e0>  " 


-r\ae  , 

j  'o,oy 


"V  - 


WJ 


which  shows  the  dependence  of  the  tracking  errors  on  the  gradients  of  the  point- 
spread  functions  at  the  origin.  In  section  6  we  shall  evaluate  these  expressions 
for  particular  array  configurations  but  first  we  give  a  short  description  of  a 
numerical  routine  to  find  targets  in  a  scene  which  is  used  together  with  the 
analytic  results  in  this  report  to  illustrate  some  array  properties. 


5  NUMERICAL  SIMULATION 


A  numerical  algorithm  has  been  written  which  takes  a  complex  image  vector 
and  searches  the  scene  for  point  sources.  The  solution  for  the  amplitudes, 
positions  and  phases  of  the  sources  is  the  one  which  minimises  the  error  between 
the  image  and  the  image  of  the  solution  (equation  (1)).  This  is  obtained  by 
using  a  nonlinear  least  squares  algorithm  (the  Levenberg-Marquardt  method)  with 
a  starting  value  provided  by  first  performing  a  coarse  search  in  position  for 
the  sources.  In  order  to  obtain  the  rms  error  in  position  for  a  single  target, 
repeated  solutions  are  made  using  a  different  noisy  image  each  time.  The  noise 
is  generated  by  using  a  pseudo-random  number  generator  to  produce  a  Gaussian 
distribution  and  hence  a  noise  vector  which  is  added  to  the  image  without  noise 
after  each  solution.  The  rms  error  is  estimated  on  the  basis  of  500  solutions. 
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EXAMPLES 


a.  Linear  Array 

Consider  an  array  of  (2K+1)  elements  positioned  at  6 
each  with  point-spread  function 

_  sin(iir-e) 

pi(6) 


0,  ±  ir,  . . .  t  ±  £ir 


(16) 


so  that  the  elements  are  at  the  Nyquist  separation  (Figure  1).  For  a 
linear  array,  the  mean  square  error  in  position  for  high  signal-to-noise 
ratio  is  given  by 


<ee>  - 


o2/2 


iKKIg] 

1 

3^ 


(17) 


which  may  be  evaluated  for  a  target  at  the  array  centre,  using  (16)  to  give 


>  ■ 


(18) 


Thus,  the  rms  error  is  inversely  proportional  to  the  square  root  of  the 
signal-to-noise  ratio  (defined,  according  to  section  3,  as  p^(0 )/cj2  «  o"2) 
and  decreases  as  the  number  of  elements  increases. 


The  variation  of  rms  error  in  position  with  target  position  can  be 
determined  by  evaluating  (17).  Figures  2a-c  plot  /(<e0>/o^)  as  a  function 
of  target  position  for  3,  5  and  7  element  arrays.  Only  the  region  6  >  0 
is  plotted  as  the  function  is  symmetrical  about  the  array  centre.  The 
positions  of  the  elements  are  at  0,  Kir.  The  figures  show  that 

the  error  is  approximately  constant  while  the  target  is  within  the 
bounds  of  the  array  but  increases  sharply  outside  the  array. 


The  dotted  lines  in  Figures  2a-c  illustrate  the  situation  where  the 
two  elements  nearest  to  the  position  of  the  target  are  used  to  estimate 
target  position,  rather  than  all  the  elements.  There  is  a  degradation  in 
tracking  accuracy,  particularly  at  the  points  where  the  target  position 
moves  from  one  pair  of  points  to  another,  though  this  is  not  very  signifi¬ 
cant.  In  practice,  in  the  presence  of  noise,  it  may  not  be  known  which 
are  the  two  array  elements  nearest  to  the  target  position,  particularly 
when  the  target  is  near  a  boundary  of  the  region.  Therefore,  it  may  be 
expected  that  the  error  is  worse  than  that  shown.  For  a  seven  element 
array  (Figure  2c),  the  ratio  of  the  rms  error  at  the  origin  using  pairs 
of  elements  to  that  using  all  the  elements  is  1.65.  Thus,  tracking 
using  pairs  of  elements  could  require  up  to  an  additional  4.35  dB  of 
signal~to~noise  ratio  to  give  equal  performance  to  that  of  the  whole  array. 

(NB.  The  expression  given  by  equation  (17)  is  valid  provided  o^/I  p?(0o)«l 
However,  this  approximation  breaks  down  at  points  outside  the  array  J 
where  Z  p?(0  )  “  0,  ie  at  multiples  of  n  for  the  sin(x)/x  point-spread 
function  nsed  here.) 


Normalised  rms  error  V«*#V**)  for  high  signal-to  noise  ratio, 
plotted  as  a  function  of  9  t  the  distance  from  the  centre  of  the 
array  for  a)  3*  to)  5#  c)  7  element  linear  array. 

The  crosses  on  the  B  axis  indicate  positions  of  array  elements.  The 
solid  curve  is  the  error  obtained  by  processing  all  array  outputs;  the 
dotted  line  that  by  processing  the  two  nearest  the  target  position. 


b.  Square  Lattice 


♦  ♦♦♦♦♦♦♦♦ 


♦  ♦♦♦♦♦♦♦♦ 

FIG.  3 


For  an  array  of  receivers  with  circularly  symmetric  point-apread 
functions  given  by 

Pi  C© .  4>)  *  J^Ort/r 

where  r2  -  02  +  <|>2,  positioned  on  a  square  lattice  at  points  (0,<|>)  ■  Gnir.nir) 
-l  <  m,n  <  i,  (see  Figure  4)  the  rms  error  in  position  has  been  evaluated 
as  a  function  of  array  size  and  target  position.  This  is  shown  in  Figures 
4a-c  where  /(<e2>/o2)  is  plotted  as  a  function  of  6  for  various  value  of  q> 
for  3x3,  5x5  and  7x7  arrays  centred  on  the  origin.  For  a  source 
which  lies  within  the  boundaries  of  the  array  the  error  varies  little 
across  the  array  but  increases  outside.  Figure  5  plots  the  error  as  a 
function  of  signal-to-noise  ratio  for  a  7  x  7  and  3x3  array  for  a  target 
at  the  origin  and  compares  it  with  the  error  produced  by  using  the  four 
nearest  elements  (taken  to  be  at  (0,0),  (0,ir),  (ir,0)  and  (tt, tr) >  to  estimate 
position.  Figures  4  and  5  show  that  for  a  2  x  2  sub-array  to  maintain  a 
constant  tracking  error  over  the  array  equivalent  to  the  larger  7x7  array 
then  up  to  4.5  dB  extra  signal-to-noise  ratio  is  required.  For  low  signal- 
to-noise  ratios,  it  is  expected  that  the  differences  in  performance 
between  the  arrays  would  be  small. 

c.  Hexagons 


F I G.  6  a 


FIG.  6b 


It  has  been  shown  elsewhere  (Webb  and  Hearn,  1984;  Petersen  and 
Middleton,  1962)  that  the  optimum  sampling  lattice  for  a  two  dimensional 
isotropic  function,  such  as  J^(r)/r,  is  based  on  a  lattice  of  equilateral 
triangles  (Figure  6a).  This  allows  the  image  plane  to  be  sampled  with  a 
lower  density  of  elements.  Results  are  presented  here  for  an  hexagonal 
array  of  seven  elements  (Figure  6b).  The  variation  of  the  rms  errors 
with  target  position  are  plotted  in  Figures  7a-c  as  a  function  of  6  for 
various  values  of  6.  As  might  be  expected,  there  are  differences  between 
the  figures  for  <e§>  and  <e?>  due  to  asynmetry.  Also  in  Figure  7  are 
plotted  results  of  the  numerical  simulation  evaluated  for  a  signal-to- 
noise  ratio  of  2.5  x  10.  Again  there  is  very  good  agreement  between  the 
simulation  and  the  theory  for  high  signal- io-noise  ratio. 


FIG.5  RMS  ERROR  IN  POSITION  AS  A  FUNCTION  OF  SIGNAL-TO  -  NOISE  RATIO  FOR 
THE  TARGET  ON  THE  ORIGIN  OF  A  7x7  ARRAY.  FOR  COMPARISON,  THE  ERROR  OBTAINED 
USING  A  3x3  AND  A  2x2  SUB -ARRAY  IS  ALSO  SHOWN~ 


riOORS  7 

Noroalisod  errors:  a)  /(Ctj >/*■*)  b)  ^  e)  /(<£?>/a'J 

for  s  sorts  tlooost  haxafonal  array  as  a  function  of  9  for  rarious  raluts  of 


7  DISCUSSION 


UNLIMITED 


The  tracking  error  for  a  single  target  as  viewed  by  an  array  of  receivers 
has  been  determined  analytically  as  a  function  of  position  for  high  s£gnal-to- 
noise  ratios.  These  results  have  been  supplemented  by  a  Monte-Carlo  simulation 
which  uses  a  numerical  algorithm  (Levenber g-Marquard t  method)  to  solve  repeatedly 
the  set  of  nonlinear  simultaneous  equations  and  estimate  the  rms  error  in 
position.  The  results  show  that  using  more  than  the  two  receiver  outputs 
allowed  by  monopulse  then  the  tracking  error  may  be  maintained  at  an  approxi¬ 
mately  constant  level  across  the  plane  of  the  array,  but  the  error  does  not 
reduce  significantly  as  more  elements  are  processed.  The  increase  in  error 
due  to  processing  two  elements  only  corresponds  to,  at  most,  a  4.5  dB  reduction 
in  signal-to-noise  ratio. 

The  numerical  method  has  been  developed  to  search  a  scene  for  many  targets 
and  Figures  8a,  b  show  20  separate  solutions  for  a  single  target  and  two-target 
scenes  for  a  seven  element  hexagonal  array  (see  section  5c) .  The  signal-to- 
noise  ratio  is  p2(0,0 )/o2,  which  is  1/4cj2  for  the  spherically  symmetric  beam- 
shape  J^(r)/r  (2.5  x  102  for  a  noise  power  of  10“3). 

Further  work  in  this  area  will  concentrate  on  analysing  the  outputs 
produced  by  more  complex  scenes  and  the  ability  of  arrays  to  resolve  unwanted 
sources.  In  addition,  optimising  the  point-spread  function  and  element  spacing 
for  single  target  tracking  will  be  considered. 

8  CONCLUSIONS 

Processing  the  outputs  from  an  array  of  receivers  does  have  definite 
advantages  over  processing  pairs  of  elements  only.  It  allows  a  target  to  be 
tracked  smoothly  across  the  plane  of  the  array,  reducing  the  error  in  the 
neighbourhood  of  elements.  An  array  based  on  a  structure  of  equilateral 
triangles  gives  similar  results  to  one  based  on  a  square  lattice  but  covers  a 
wider  field  of  view  for  a  given  number  of  receiver  elements. 
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